Strength Reduction in Electrical and Elastic Networks 
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Particular aspects of problems ranging from dielectric breakdown to metal insulator transition can 
be studied using electrical or elastic networks. We present an expression for the mean breakdown 
strength of such networks. First, we introduce a method to evaluate the redistribution of current 
due to the removal of a finite number of elements from a hyper-cubic network of conductances. It 
is used to determine the reduction of breakdown strength due to a fracture of size k. Numerical 
analysis is used to show that the analogous reduction due to random removal of elements from 
electrical and elastic networks follow a similar form. One possible application, namely the use of 
bone density as a diagnostic tools for osteoporosis, is discussed. 
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I. INTRODUCTION 

Networks of electrical or elastic elements are used as 
models to study phenomena observed in disordered sys- 
tems. They include dielectric breakdown pi 0, 9, metal- 
insulator transitions Q , brittle fracture in disordered 
solids |], [| fl], ||] , and strength of trabecular bone || . In 
particular, such models have provided insights on the 
critical phenomena [|l0| , scale- invariant disorder O, [To) , 

In this paper, 
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Breakdown processes depend on the size of the longest 



and size-dependence of the system 
we will study variations in strength of a class of such 
systems. 

:akdown pro c 

defect |l2|, [l^, ]l4| and such extreme value problems are 
difficult to analyze |n^. However, their statistical proper- 
ties are expected to be amenable to analysis; this is seen 
to be the case for the class of problems we study here. 

We first consider fused-conducting networks, where the 
breakdown currents of each element is assumed to be pro- 
portional to its conductance; i.e., breakdown of a conduc- 
tor occurs when the potential difference across it reaches 
a fixed value. In the simplest case, conductances along 
each axis of a hyper-cube are set equal. In Section II, we 
calculate how an external current introduced at a node 
and removed from an adjacent node is distributed on the 
network |l5|, [Tq ]. The Green's Function derived for this 
case can be used to calculate the current distribution 
due to an external field on a network from which a finite 
number of conductances are removed (Section III). This 
calculation is used to determine how the breaking current 
on a network of fused conductances is reduced due to a 
fracture of size n. The remaining sections of the paper 
involve numerical analysis to show how this relationship 
is modified in a variety of other electrical and elastic net- 
works. Section III also provides a discussion of the effects 
of random removal of elements from a two dimensional 
network. 

In Section IV, we present analogous results for cubic 
networks with randomly chosen conductances. The rela- 
tionship between the strength of a network and the frac- 
tion of bonds removed is shown to have a form similar to 



relationships found in section III. Section V shows that 
the strength of disordered elastic networks which include 
elastic and bond-bending energies also follow the same 
form. 

The motivation for this study is to develop a method 
to relate the mean strength of a bone to its mass [|l7], [l8| . 
A solution of this problem can provide a non-invasive 
method to diagnose bone strength and can aid in the 
management of osteoporosis. These issues are addressed 
in the concluding section. 



II. GREEN'S FUNCTIONS FOR HYPER-CUBIC 
NETWORK 

Consider an infinite, d-dimensional hyper-cubic net- 
work where all conductances along the m th direction 
(u m ) are assumed to be equal to a m (m = 1, 2, d) . In 
this section, we will calculate the current distribution on 
this network due to a unit external current introduced 
at the origin and removed from an adjacent node, say 
a = — fli . This result was calculated using a Green's 
function method by Kirkpatric JTT| , |l5|] , and is reproduced 
here for completeness. 

Denote the potential at the node n = (m, ri2, rid) , 
by V(n) and the current on the conductor joining nodes 
n and n + u m by J m (n) ; then, 



J m (n) = a m [V(n) - V(n + u m )] 



(1) 



The potentials V(n) can be solved using the Kirchhoff's 
rules at each node; i.e., 



Ed 
m=l a * 



{[V(n) 
[V(n) 



y(n + u m )] 
F(n-u m )]} 



(6» 



<5n,a) 



(2) 



where the right side is the externally applied current. 
Eqn. (2) is easily solved by using the Fourier transform 
Vk = J2n e ~ mk ^n ; which satisfies 



2 a m [1 - cos(u m • k)] V-(k) = (1 - e 



ia k\ 
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Thus 

where k m = u m 
V(n) 



(1- 
Em=l 2cr « 
k. Hence 

1/2 



and 
Jm(n) 



cr m /2 r7T 
_ (27r) d J-tt 



<ik e 



e -iak^ 

[1 - cos(fc m )] 



dk e mk V(k), 



ink 



(3) 



(4) 



(5) 



El, = l<Tm/[l-COs(fc m ,)] 



Next we specialize to networks in two dimensions with 
equal conductances (<x m = 1), and define a(n x ,n y ) = 
J y (n x ,n y ) and (3{n x ,n y ) = J x (n x ,n y ) , see Figure 1. 
The symmetries of the network imply that 



a{n x ,n y ) = a(-n x ,n y ) = a(n x ,-n y ), 
f3(n x ,n y ) = -(3(-n x ,n y ) = -f3(n x ,-n y ) . 



(6) 



Combining Eqns. (5) and (6) give 

2 f 
~ ^ JO 



a(n x , n y ) 
/3(n x ,n y ) 



7? f( 



cos(n x k x ) cos(n.y k^) ■ 9 / kv \ 11 

— — -- sin {-$-) dk, 

1 I 

y* 2 

ky) 

x sin(^)sin(^)dk. 



2 — cos(k x )~- cos(k y ) 
7T sin [(n x + ^)kx\ sin [(n y -\-^)k y ] 



(7) 



2— cos(fe x 



-cos(k y ) 



Values of a's and /3's for several pairs of (n x ,n y ) are 
presented in tables I and II, and will be used for compu- 
tation in Section III. Notice that a(0,0) = 1/2, as can 
be easily confirmed using arguments based on symmetry 
and superposition pa] . 
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FIG. 1: Scheme to compute the current distribution due to 
a unit external current introduced at the origin and removed 
from an adjacent node. 



TABLE II: Values of j3 for the first pairs of (n x , n y ) . 



n v 





1 


2 


3 


4 


5 





0.18169 


0.04507 


0.01315 


0.00469 


0.00205 


0.00106 


1 


0.04507 


0.03052 


0.01596 


0.00826 


0.00451 


0.00264 


2 


0.01315 


0.01596 


0.01192 


0.00789 


0.00510 


0.00334 


3 


0.00469 


0.00826 


0.00789 


0.00626 


0.00464 


0.00337 


4 


0.00205 


0.00451 


0.00510 


0.00464 


0.00384 


0.00304 


5 


0.00106 


0.00264 


0.00334 


0.00337 


0.00304 


0.00259 


6 


0.00062 


0.00165 


0.00225 


0.00244 


0.00236 


0.00214 



TABLE I: Values of several a(n x , n y )'s. In particular 
a(0,0) = 1/2, and q(1,0) = 1/2(4/tt - 1). 



7ly_ 





1 


2 


3 


4 


5 


o 


0.50000 


0.13662 


0.04648 


0.02019 


0.01080 


0.00670 


1 


0.13662 


0.00001 


0.01455 


0.01174 


0.00816 


0.00571 


2 


0.04648 


0.01455 


0.00000 


0.00404 


0.00441 


0.00383 


3 


0.02019 


0.01174 


0.00404 


0.00000 


0.00162 


0.00208 


4 


0.01080 


0.00816 


0.00441 


0.00162 


0.00000 


0.00080 


5 


0.00670 


0.00571 


0.00383 


0.00208 


0.00080 


0.00000 


6 


0.00457 


0.00413 


0.00314 


0.00205 


0.00113 


0.00045 



the y— direction. The currents on this configuration are 
Jx(nx, n y) = , J y {n x ,n y ) — 1, see Figure 2(a). In this 
section , we introduce a method to calculate the current 
distribution on the network when a finite number of con- 
ductances are removed. As an application we determine 
the enhancement of current on the edges of a "fracture" 
formed by removing n consecutive conductances. 



A. Removal of a Single Conductor 



The a's and /3's can also be computed perturbatively 
by expanding the denominator as a series in terms of 
(cos k x + cos k y ) /2. Expansion to fourth order gives val- 
ues for a(n x ,n y ) and f3(n x ,n y ) that are accurate to 
within two percent. 

III. RESISTOR NETWORKS IN TWO 
DIMENSIONS 

Consider next, a two dimensional isotropic (a x = a y = 
1) network with an externally applied unit electric field in 



First, we enumerate changes in the current distribu- 
tion due to the removal of the conductor joining the ori- 
gin to (0, —1), Figure 2(b). The excess current 5 on any 
remaining element can be considered to be the redistri- 
bution of the current that originally passed through the 
link (0,0) — (0, —1), see Figure 2(c). 5 can be evaluated 
using the values of a(m, n) and [3(m, n) by the following 
procedure. Consider the complete network, with and ex- 
ternal current (1 + J) introduced at (0,0) and removed 
from (0, —1), as shown in Figure 2(d). (1 + J) is cho- 
sen so that a current J passes through (0,0) — (0, — 1), 
and the remainder (= 1) passes through the rest of the 
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FIG. 2: Scheme to compute the current distribution when 
a single conductance is removed. In sequence: (a) an isotropic 
network with a current flow along the y— direction. In (b) 
a conductance is removed , and 8 is the excess of current in 
any remaining element. Next in (c) , a unit external current 
is introduced in order to compute 8. Finally (d) shows the 
configuration from which 8 can be evaluated. 



network; i.e., this second part is the solution to the prob- 
lem illustrated in Figure 2(c). But from the discussion 
in section II, the current passing through (0, 0) — (0, —1) 
is q(0, 0) x (1 + J). Hence, 



J = a(0,0) x (1 + J) 



(8) 



Since a(0,0) = 1/2, we find that J = 1 and that the 
changes of current in the network are given by, 



AJ y (n x ,n y ) = (1 + J)a(n x ,n y ) = 2a(n x ,n y ) 
AJ x (n x ,n y ) = (1 + J)f3(n x ,n y ) = 2 (3(n x , n y ) . 



(9) 



In particular , the current on conductances adjacent to 
that removed in Figure 2(b) is 



l + 2a(l,0) = 4/tt, 



as was given by Duxbury et. al. 

As a second illustration, we calculate the currents on 
the network shown in Figure 3(a). As in the first exam- 
ple, the problem can be solved using a complete network 
on which external currents shown in Figure 3(b) are ap- 
plied. Notice that prior to removal of the conductances 
there are no current in AB and hence the current J\ that 
is removed at B is required to pass entirely through AB. 
The currents J\ , J 2 , J 3 can thus be obtained from 



Ji = a(0,0)Ji + 0(0,0) (1 + J») + /3(1,0)(1 + J 3 ), 

J 2 = /3(0,0)Ji + a(0,0)(l + J 2 ) + a(l,0)(l + J 3 ) , (10) 
J 3 = /9(-l,0)Ji + a(-l,0)(l + J 2 ) + 0(0,0) (1 + J 3 ). 



Notice that the roles of a's and /3's are reversed when 
computing the distributions due to the external current 
J 3 . Solving these equations give Ji = 1.50, J2 = 2.38, 
and J 3 = 2.06 . Some of the currents on this network 
are given in Figure 3(a), and were confirmed through 
numerical integrations. 



Finally we consider the distribution of current due 
to a "fracture" of k consecutive cr^'s lying along the 
a;— direction, and compute the largest currents (i.e., those 
on the conductances forming the edges) on the network. 
Assume that the conductances between (i,0) — (t, — 1) 
where t — 1, 2, n arc removed. This problem is equiv- 
alent to having a complete network with currents ( Jt + 1) 
introduced at node (t,0) and removed at (t, — 1) such 
that currents Jt pass through links (t, 0) — (i, — 1). These 
external currents can be obtained by solving the set of 



equations 

K 

J t = 5^a(B-*,0)(J. + l), (11) 

s=l 

and the current on the conductance (0,0) — (0,-1) is 
given by 

K 

J max (K) = l + ^a(s,0)(J s + l). (12) 

s=l 

We have solved Eqns. (11) for values of k between 1 
and 250, and evaluated the maximum current Jmax{^) on 
the network using Eqn. (12). Their values as a function 
of k are shown by circles in Figure 4. To a very good 
approximation, J max is given by 

J max (K) w 1 + 01 k 1/4 + a 2 k 1/2 , (13) 



4 



with ai ~ —0.582 , and a-2 ~ 0.894 ; this expression is 
shown by the solid line. For k » 1, this approaches the 
form Jmax{n) = Jo (l + a k 1 / 2 ) , proposed in Ref . [pd| . 
Expressions (11), (12) and (13) were confirmed numeri- 
cally in a finite network of size L x x L y (= 100 x 100) cen- 
tered at the origin. For k > 45 effects of the boundaries 
become significant and values of J max near the bound- 
aries are larger than those given by Eqn. (13). 

Next, consider a square network of fused conductances; 
i.e., if the current on an element exceeds To, then it will 
"burn" , and its conductance would become zero irre- 
versibly. The complete network will fuse when the ex- 
ternal field is increased to E = I max {0) = Iq, when the 
current on each element in the y— direction becomes Jo- 
The presence of a fracture of length k will reduce the 
current I m ax(n) per link introduced at the top edge of 
the network. To calculate I m ax{n), notice that the first 
elements on the network to burn will be those on the 
edges of the fracture. As a result k will increase, and 
the fracture will propagate in either direction, until the 
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FIG. 3: (a) An isotropic network with the three conduc- 
tances removed, and current distributions on its boundary. 
The scheme to compute 5 is shown in (b). 



FIG. 4: The largest current distribution resulting from 
a removal of n consecutive conductances lying along the 
a;— direction. Results obtained by using Eqns. (11) and (12) 
are shown by circles, while the expression (13) is shown by 
continuous line. 



network is separated into two segments. From Eqn. (13) 

^max(^) 1 



,(0) 



1 + aiK c 



a 2 K 



2d 



(14) 



with a = 1/4 . In the remainder of this section an the 
next, we will show through numerical analysis, how this 
form is modified when resistors are randomly removed 
from square and cubic networks. Further, the failure 
strength of random resistors networks and disordered 
elastic networks will be shown to exhibit a similar form, 
as long as the networks are not close to the bond perco- 
lation threshold. 

We also note that the maximum current I m ax does not 
necessarily correspond to the first fracture of a conduc- 
tance in these random cases O, 19 . However, the yield 



and ultimate currents have been shown to exhibit similar 
behavior. 



B. Random Removal of Conductances from a 
Single Layer 



Here we consider a finite network of size L x x L y 
and study the reduction of the maximum current, say 
Imax{ v ) , due to the random removal of a fraction v of 
conductances on a single layer y — . The length of a 
largest fracture k on the layer can be estimated by pl| 



L x v K (1 



1 



(15) 



5 




0.6 - 




FIG. 5: The reduction of Imaxiy) due to the random removal 
of a fraction v of conductances on a single layer y — follows 
Eqn.(17). Here we show the average and the standard error 
of a set of 100 networks for values of v < 0.70, with a ps 
0.37, oi « -0.90 and a 2 w 1.17. 



FIG. 6: The reduction of the maximum current as a func- 
tion of the fraction of conductances removed from a square 
network. We considered ensembles of 25 networks with size 
100 x 100. Expression (18) is plotted by continuous line and 
the average value of a ~ 0.55, ai « 0.09 and 02 « 0.20. 



where the probability on the left is obtained by requiring 
the need for k absent locations (i.e., v K ) bounded by two 
links (i.e., (1 - v) 2 ). Thus if (1 - v) 2 « L x 



K ~ log(V")' 
Assuming a form similar to Eqn. (14), 



c(0) 



(16) 



(17) 



Figure 5 shows results from computations on a set of 
100 networks, for values of v < 0.70 . In these integra- 
tions the boundaries at the top and bottom were kept at 
potentials Vq and 0, respectively. Had the fractures been 
independent of each other we would expect a = 1/4; the 
fact that a « 0.37 suggests a strong correlation between 



fractures [19 



follows a relationship similar to Eqn. (17) where a de- 
pends on e . 

Consider first the case e = 0; i.e. all conductances 
are equal. As the fraction v of conductances removed 
reaches the bond percolation threshold up (= 1/2), the 
network separates into multiple segments (EOl Elf , and 
hence J mol (i/o) — ► 0. Eqn. (17) fails to satisfy this con- 
dition. (Notice that for the quasi-lD case discussed in 
Section IIB, v$ = 1). We propose a modification 



1 



,(0) 



i+„ r log w \ a , r iog(jv) i 
1 + ai \i°im; +fl2 iM7)i 



(18) 



L y is the 



of Eqn. (17) for this case. Here N = L x x 
number of nodes on the network. Figure 6 shows that the 
mean values (of 25 networks of size 100 x 100) satisfies 
Eqn. (18) with a value of a s=s 0.55 . Including random- 
ness in the conductances changes the value of a, a\ and 
CI2; but de agreement with Eqn. (18) is equally good . 



Random Networks in Two Dimensions 



IV. RESISTOR NETWORKS IN THREE 
DIMENSIONS 



We extend the computations of Section IIB to a net- 
work of random fused conductors from which a fraction 
v of elements have been removed. The conductances are 
chosen randomly within a range [1 — e, 1 + e] and fusing 
current of each element is assumed to be proportional to 
its conductance (i.e., they burn when the potential differ- 
ence between the ends exceed a fixed amount). Imaxiy) 



In this section, we discuss the behavior of the break- 
ing strength of cubic networks of fused conductances. 
Green's function calculations (of Section II) can be used 
to evaluate the current distributions via a(n x ,n yi n z Ys 
and P(n x ,n y ,n z ys. Changes in the current distribution 
due to the removal of a finite number of conductances 
can be evaluated using methods outlined in section IIIA. 



6 



TABLE III: Values of a's on the plane z — , for pairs of 

(Tlx 5 ) ■ 



tly 





1 


2 


3 


4 


5 





0.33333 


0.06175 


0.01392 


0.00404 


0.00153 


0.00073 


1 


0.06175 


0.02323 


0.00793 


0.00298 


0.00130 


0.00066 


2 


0.01392 


0.00793 


0.00379 


0.00184 


0.00095 


0.00054 


3 


0.00404 


0.00298 


0.00184 


0.00109 


0.00066 


0.00041 


4 


0.00153 


0.00130 


0.00095 


0.00066 


0.00045 


0.00031 


5 


0.00073 


0.00066 


0.00054 


0.00041 


0.00031 


0.00023 


6 


0.00040 


0.00038 


0.00033 


0.00027 


0.00022 


0.00017 



A. "Penny-Shaped" Fracture 

We present the calculation of the current enhance- 
ment due to a "penny-shaped" fracture in the Z = 
plane jyj [l^] , centered at the origin. For a given radius 
p all conductances in the z— direction within a distance p 
of the origin in the x — y plane are removed and changes 
in the current distribution are computed. This calcula- 
tion requires evaluating 



a(n x , n y , 0) = -^r 



cos(n x k x ) cos(n y k y ) sin 2 (fc z /2) 
3 — cos^) — cos(fc y ) — cos(fc 2 )' 



r/k 



(19) 



Some values of a(n x , n y , 0) are given in TABLE III. The 
currents on the edge of a fracture can be calculated by 
methods similar to those in section III. It is found that 




the maximum current for this case is approximately, 

J max {K) « 1 + oi fc 1/8 + a 2 fc 1/4 , (20) 

with ai « —0.95 and a 2 ~ 1.07, see Figure 7. Once 
again these conclusions were checked in numerical inte- 
grations on cubic networks of conductances. 

B. Random Resistor Networks 

In these computations a fraction of conductances from 
an L x x L y x L z network was removed randomly. A 
potential Vq was imposed at the top surface of the net- 
work, and the bottom surface was grounded. The be- 
havior of the maximum current Ima X {y) for increas- 
ing v is recorded. As in square networks, I ma xiy) is 
seen to satisfy Eqn. (18), where for cubic networks 
u w 0.7508 & pn). As before, it is found that a de- 



pends on the level e of disorder. 



V. ELASTIC NETWORKS 

The elastic network model we consider is constructed 
by displacing the nodes of a hyper-cubic network ran- 
domly |J . The displacements (of less than e ) are chosen 
to be small enough that the topological structure of the 
network remains unchanged. Adjacent nodes are joined 
by struts whose elastic moduli are chosen randomly from 
a pre-determined range; i.e., k e € [&o(l — rj e ) , ko(l + r] e )} . 
Bond-bending terms which depend on the variation of 
the angles between adjacent bonds on the network are in- 
cluded. The linear response coefficients are chosen from 
within a range Kb G [«o(l — ? 7b),Ko(l + ?/;,)] . These 
networks are degraded by removing a fraction v of the 
links randomly, and the remaining elastic elements are 
assumed to retain their moduli. Thus the total energy of 
the network is 



U 



elastic 
elements 



\k e {8xf+ ]T \n b (59t 



bonds 
angles 



where S x is the change of the length of an elastic element 
under a given external strain of the network, and 89 is 
the variation of a bond-angle from equilibrium. 

A strain based criterion is used to implement fracture 
of elastic elements and bond angles of the network. Any 
strut that is strained beyond a value 7 is removed from 
the network along with all bond-angles it contributes to. 
7 for the elastic elements are chosen from a Weibull dis- 
tribution whose cumulative density is p2[ 



FIG. 7: The largest current resulting from a removal of a 
"penny shaped" fracture of radius p lying on the plane z = 0. 
The calculations using Eqn. (19) are shown by circles, while 
the expression (20) is shown by the line. 



(7(7) = 1 — exp 



(21) 



A corresponding criterion is included for fracture of 
bonds. 



0.0 0.1 0.2 0.3 0.4 0.5 

V 

FIG. 8: The reduction of T m ax[v) due to the random removal 
of a fraction v of elastic elements. In continuous line is shown 
the average obeying expression (22) of a square network of size 
40 x 90. The parameters for the best fit are a m 0.97, a± « 
-0.11, a 2 « 0.07. 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

V 

FIG. 9: The same as in Figure 8 for a cubic elastic network of 
size 20 x 20 x 30. For the nonlinear fit a w 1.14, ai ~ -0.011, 
and a-2 ~ 0.015. 

VI. DISCUSSIONS AND CONCLUSIONS 



A given uniform external strain £ is imposed through 
a sequence of small adiabatic increments. If elastic ele- 
ments or bond-angles are removed from the network due 
to a fracture, equilibrium is re-calculated prior to increas- 
ing £. Nodes on the sides of the network transverse to 
the external strain are constrained to move in plane. The 
force T(£) required to sustain a strain £ is calculated by 
adding the vertical forces on the top of the layer. The 
equilibrium of a network is calculated by minimizing its 
energy using the conjugated gradient method pi| . 

For small values of the compression £, the stress-vs- 
strain relationship T(£) is linear. Numerical analysis 
shows that T(£) becomes nonlinear (at the yield point) 
following the first fracture of elastic elements || . As £ is 
increased further, the stress reaches a maximum (which, 
for a network a fraction v of whose elements have been 
removed, is denoted T max (u)) prior to failure. 

The analog of Eqn. (18) for an elastic network is 

Tmaxiy) _ 1 (22) 

T max(0) f log(iV) Y , f log(jV) ] 

Figures 8 and 9 show that this expression is satisfied for 
square and cubic disordered elastic networks. As in elas- 
tic networks, a depends on the quantities e, rj e , and rjb 
characterizing the level of disorder. It is further found 
to depend on the ratio k e / Kb of mean linear response 
coefficients of elastic and bond-bending terms. 



The motivation for the work reported here is the need 
to identify measures that can be used as non-invasive di- 
agnostic tools to estimate the strength of bone. Large 
bones consist of an outer compact shaft and an in- 
ner porous region (trabecular architecture) whose struc- 
ture is reminiscent of a network of disordered cubic net- 
works |J . Due to loss of fracture toughness jL4j of the 
outer shaft, the trabecular architecture becomes the prin- 
cipal load carrier in older adults. Reduction in sex- 
steroid hormones (estrogen and testosterone) effect bone 
turnover, leading to degradation of the porous bone. Be- 
yond this stage, it is important to monitor the loss of 
bone strength at regular intervals, preferably with non- 
invasive diagnosis tools. 

The principal mechanism for loss of trabecular mass 
is through removal of individual struts, due to minor 
trauma. Fracture toughness and thickness of the remain- 
ing trabecular bone do no change appreciably. Mechan- 
ical studies on ex-vivo bone samples have shown that 
trabecular networks from patients with a broad range of 
ages fracture at a fixed level of strain [^3) , even though 
the corresponding fracture stresses exhibit large varia- 
tions. These observations motivated the fracture criteria 
used in our models. 

Eqn. (18) is the main result of the paper. Its form is 
justified by the calculations on a network of fused con- 
ductors presented in Section IIIA. The numerical anal- 
ysis of random electrical networks and disordered elas- 
tic networks are used to confirm that the strength bears 
the same relationship to the level of elements removed 
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from the network. The index a corresponding to ran- 
dom removal of elements is larger that in the case of a 
single fracture (=1/4), indicating that correlations be- 
tween fractures is important to determine the strength 
of a network [jl9| . 

It is important to test if samples of trabecular bone sat- 
isfy Eqn. (18). Previous studies on bone samples have 
suggested a nonlinear relationship between the breaking 
strength (T) and the effective density (p) of trabecular 
bone. It was found heuristically that T » p" 13 , where 



(3 w 2.6. If it is assumed that the density of bone ma- 
terial remains intact and that the fractional reduction of 
effective density is proportional to is, then the relation- 
ship Eqn. (18) can be tested using this data. Such a 
comparison will be done in future work. 
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